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Abstract 


Leighton Chajnantor Telescope (LCT), i.e., the former Caltech Submillimeter Observatory telescope, will be 
refurbished at the new site in Chajnantor Plateau, Chile in 2023. The environment of LCT will change significantly 
after its relocation, and the telescope will be exposed to large wind disturbances directly because its enclosure will 
be completely open during observation. The wind disturbance is expected to be a challenge for LCT’s pointing 
control since the existing control method cannot reject this disturbance very well. Therefore, it is very necessary to 
develop a new pointing control method with good capability of disturbance rejection. In this research, a disturbance 
observer—based composite position controller (DOB-CPC) is designed, in which an Ha feedback controller is 
employed to compress the disturbance, and a feedforward linear quadratic regulator is employed to compensate the 
disturbance precisely based on the estimated disturbance signal. Moreover, a controller switching policy is 
adopted, which applies the proportional controller to the transient process to achieve a quick response and applies 
the DOB-CPC to the steady state to achieve a small position error. Numerical experiments are conducted to verify 
the good performance of the proposed pointing controller (i.e., DOB-CPC) for rejecting the disturbance acting 
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on LCT. 
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1. Introduction 


The Caltech Submillimeter Observatory (CSO) telescope 
was decommissioned in 2015 and will be moved to its new 
site at Chajnantor Plateau, Chile from its old site at 
Maunakea, Hawaii in 2023. After the completion of the 
relocation, the CSO telescope will be refurbished at the new 
site and renamed as the Leighton Chajnantor Telescope 
(LCT) (Vial et al. 2020; Chen et al. 2022). However, the 
relocation and refurbishment of LCT will face with new 
challenges for its pointing control system, because it is 
expected to resist strong wind disturbance at the new site, 
which has an average wind speed of 10 ms ', a maximum 
instantaneous wind speed of 30 ms‘, and more frequent 
wind speed fluctuation (ALMA 2010). In addition, LCT will 
work with its enclosure open completely, because LCT is 
required to observe a larger area of the sky at the new site to 
achieve the new scientific goals (Caltech et al. 2016). 

To achieve a high pointing accuracy at the new site under 
strong wind load, LCT’s pointing control system needs to 
realize a pointing control with high precision, i.e., to ensure 
the azimuth angle and the elevation angle of LCT’s antenna 
be within a sufficiently small deviation from position of the 
target source in the sky. Since LCT is a submillimeter radio 
telescope with a 10.4 m-diameter antenna and receives the 
radiation with the wavelengths from 2 mm to 350 um, a 


pointing error below 3” (i.e., 0°000 833) is required to 
achieve its scientific goals,” which requires that the root mean 
square (rms) of the position error of LCT’s antenna must be 
less than 3”. However, this pointing accuracy is very difficult 
to achieve when LCT is working with a strong wind load. By 
testing the existing proportional feedback controller of LCT 
on its simulation model (Chen & Wang 2022), the rms of the 
azimuth angle error of LCT’s antenna is greater than 3” when 
the wind with a speed above 10 ms’! is introduced. 
Obviously, the performance of the existing feedback 
controller cannot satisfy the requirement on the pointing 
accuracy when the telescope is fully exposed to high-speed 
wind. Thus, for the purpose of keeping the pointing accuracy 
of LCT at the new site, it is necessary to improve the position 
controller of LCT. 

Few work has been done to employ advanced control 
methods to the pointing control systems of astronomical 
telescopes. In practice, traditional proportional-integral (PI) 
controller is the most popular one applied to the pointing 
control systems of telescopes due to its convenience for 
implementing (Ravensbergen 1994; Burns 1995; Gawronski 
2001; Ranka et al. 2016), although it is sensitive to wind 
disturbance. In the past decades, to make the position controller 
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of radio telescopes more accurate and robust under external 
disturbance, many scientists and engineers have paid a lot of 
efforts on improving its performance. On one hand, some 
model-free control methods are developed since they have no 
strict requirements on the accuracy of the model of the 
telescopes’ pointing control systems. For example, the active 
disturbance rejection controller (ADRC) is designed based on 
the classical feedback control framework. With the help of 
tracking differentiator, feedback control law, and the extended 
state observer (ESO), the ADRC controller can resist the 
internal uncertainty (e.g., the nonlinearity of the system) and 
the external disturbance (e.g., the wind disturbance) (Qiu et al. 
2014; Ranka et al. 2015; Guo et al. 2016; Feng & Guo 2017; 
Deng et al. 2018; Li et al. 2019). However, as the core of 
ADRC controller, the ESO can only estimate rough disturbance 
signal and is very sensitive to parameter tuning since it is not 
designed based on the system’s model. On the other hand, since 
more methods are developed to identify the models of the 
telescopes’ pointing control systems, it is feasible to employ 
model-based control methods for improving the robustness. 
The Linear-Quadratic-Gaussian (LQG) controller is a model- 
based optimal control method that is widely applied to many 
modern radio telescopes, e.g., the Deep Space Network (DSN) 
(Gawronski et al. 1994; Maneri & Gawronski 2000), the Green 
Bank Telescope (GBT) (Gawronski & Parvin 1995), the Large 
Millimeter Telescope (LMT) (Gawronski & Souccar 2005), the 
Telescope Nazionale Galileo (TNG) telescope (Schipani et al. 
2020), and the 110m QiTai radio Telescope (QTT) (Li et al. 
2017). The introduction of the LQG controller effectively 
restrains the vibration of these radio telescopes such that the 
error of the position control of their pointing control systems 
can be reduced under wind disturbance compared to the 
position error under PI controllers. As the H norm limits the 
performance of the LQG controller, the Ha controller is 
proposed by replacing the Hz norm with the H,, norm to 
further improve the performance of the optimal controller 
(Gawronski 2001). 

Although many model-free or model-based control methods 
have been applied to telescopes’ pointing control systems, the 
effect of disturbance can hardly be eliminated since it is not 
accurately estimated. An H,, controller was designed by Chen 
& Wang (2022) for LCT’s pointing control system to improve 
its robustness for resisting wind disturbance. From the 
simulation results, although the H,, controller can significantly 
compress the position error of LCT’s antenna, it still cannot 
eliminate the wind disturbance precisely. Hence, in this paper, 
to eliminate the wind disturbance precisely, we design a 
disturbance observer (DOB) to realize an accurate estimation of 
wind disturbance, based on which we propose a new 
disturbance observer-based controller (DOBC) for LCT to 
acquire better performance of disturbance rejection. 

DOBC offers a feasible method for reducing the negative 
influence from external disturbance and internal uncertainty by 
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utilizing the estimated disturbance from the DOB and feed it 
forward to the control signal of an electromechanical system 
(Choi et al. 2003; Chen et al. 2016). DOBC has been applied to 
many real-world systems, e.g., permanent magnet synchronous 
motors (PMSM) (Dai et al. 2021), missile seekers (Sadhu & 
Ghoshal 2011), unmanned aerial vehicles (UAV) (Huang & 
Chen 2022; Tripathi et al. 2022), overhead cranes (Wu et al. 
2020), ship-mounted tower cranes (Qian & Fang 2019), cars 
with an active suspension system (Pan et al. 2016), and dual- 
flexible manipulators with a telescopic arm (Shang et al. 2022). 
For generating a compensation control signal based on the 
estimated disturbance, most research on designing DOBC 
(Sadhu & Ghoshal 2011; Pan et al. 2016; Qian & Fang 2019; 
Wu et al. 2020; Huang & Chen 2022; Shang et al. 2022; 
Tripathi et al. 2022) simply assume that the control signal (i.e., 
the output of the controller instead of the actuator) and the 
disturbance signal are of the same type (e.g., both are torques or 
both are forces). Based on this assumption, the estimated 
disturbance signal (i.e., the output of the DOB) is directly 
utilized as a compensation control signal (by taking it negative 
value) and is simply combined with feedback control signal. 
For example, Pan et al. (2016), Qian & Fang (2019), and Wu 
et al. (2020) ideally consider the feedback control signal to be 
force, and the disturbance acting on the controlled object is also 
considered as force so that the estimated disturbance signal can 
be directly utilized as the compensation control signal. 
However, in realistic control systems, the output of the 
controller is an electric signal (e.g., voltage) rather than force 
or torque. The output of the controller (i.e., the control signal) 
is sent to the actuator and the actuator transforms the control 
signal to force or torque, which directly control the object. 
Therefore, to compensate the disturbance, the output of the 
DOB, which represents some physical quantity (e.g., force or 
torque), must be transformed to a compensation control signal 
that will be sent to the actuator and transformed to the physical 
quantity (e.g., force or torque) acted on the controlled object. 
Some researchers designed a compensation gain (Chen et al. 
2016) or a nonlinear continuous function (Dai et al. 2021) to 
transform the estimated disturbance signal to a compensation 
control signal. In this paper, we construct an optimal tracking 
control model with a linear quadratic regulator to generate a 
precise compensation control signal that can make the actuator 
output proper physical quantity relatively close to the estimated 
disturbance signal such that the disturbance can be 
compensated. 

The rest of this paper is organized as follows. In Section 2, 
the model of LCT’s pointing control system is constructed with 
the vibration modals of the antenna, which are the first-order 
and the second-order vibration modals of the antenna. Section 3 
presents the design of the DOBC based on the model of LCT’s 
pointing control system. The results and analysis of the 
simulation are given in Section 4 to verify the performance of 
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the DOBC. In Section 5, we conclude this paper and give the 
outlook for future research. 


2. Model of LCT’s Pointing Control System with 
Structural Vibration Modals 


An accurate model of LCT’s pointing control system is very 
important to the design of the position controller. LCT’s 
pointing control system model was constructed in our previous 
work (Chen & Wang 2022), which consists of the position 
controller, the DC motor, the reducer box, and the antenna. 
There are two major vibration modals of the LCT antenna, 
which are the low-frequency flexible modal (i.e., the first-order 
vibration modal of the antenna), and the main structural 
vibration modal of the main reflector (i.e., the second-order 
vibration modal of the antenna). In our previous work, the first- 
order vibration modal of the antenna is already considered 
since it is the low-frequency flexible modal resulting from the 
flexible connections between the DC motor and the reducer and 
between the reducer and the antenna when the antenna is 
regarded as a rigid body. Thus, the dynamic equation of the 
first-order vibration modal can be formulated as shown in 
Equation (1). For the clarity of description, in this paper, we 
only consider the position control in azimuth. The method 
developed in this paper can also be applied to the position 
control in elevation. 


‘3 ; 1 
J,9, + D101 = Kroc Troc = Keoe( Abc a) (1) 


where J4 is the moment of inertia of the LCT antenna; D, is the 
damping of the rigid body of the LCT antenna; Kppc is the 
stiffness of the reducer, which takes the flexible connections 
among the DC motor, the reducer, and the antenna into 
account; Troc is the output torque of the reducer; @pcy is the 
rotation angle of the DC motor; and 6, is the output azimuth 
angle of the rigid body of the LCT antenna. Then, Equation (1) 
can be reformulated as 


Js; + DO; + Krochi = 


K 
y ODM (2) 


Applying the Laplace transform to both sides of Equation (2) 
yields 


s? Ja 0i + sDı0ı + Krpocôı = 


K; 
y ODM (3) 


Therefore, the transfer function of the first-order modal is 


O Kroc/N 


G(s) = z 
Opcem Jas? + Dis + Kroc 


(4) 
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Naturally, we can obtain the angular velocity w, of the first- 
order modal as 


_ | Kroc 
w= |" 
Ja 
_ {6.99 x 10’N- m- rad"! 
160560 kg - m? 


æ% 20.8651 rad-s-! (5) 


So, the natural frequency of the first-order modal is 
fi ~ 3.32 Hz, which is also verified by the modal analysis of 
LCT using Ansys 2022. 

Although flexible connections among the DC motor, the 
reducer, and the antenna are considered in the first-order modal, 
the antenna is still assumed to be a rigid body at present 
although LCT’s antenna in the real world is a flexible body 
(which may lead to the deformations of the main reflector under 
the wind and gravity). Therefore, to better describe the 
flexibility of LCT’s antenna, we introduce the structural 
vibration modal (which is related to the vibration of the main 
reflector) as the second-order vibration modal of the antenna to 
extend the pre-established antenna’s rigid body model. By 
analyzing the vibration modals of LCT’s antenna with Ansys 
2022, we find that the natural frequency of the main structural 
vibration modal of the antenna is 13.004 Hz. By introducing 
the main structural vibration modal, the extended LCT antenna 
model can be formulated as follows: 


Jaðı + D1 0, = T; 
Jaĝı + D202 + Ko. = T 
04 = 01 + O2 (7) 


where J, is the moment of inertia of LCT’s antenna in azimuth; 
D,, Dz are the damping of the rigid body and the main 
structural vibration modal respectively; K2 is the stiffness of the 
main structural vibration modal; 04, 02 are the output azimuth 
angles of the rigid body and the main structural vibration modal 
respectively; 04 is the azimuth angle of LCT’s antenna (in 
rad - s7"); and T; is the load torque of LCT’s antenna. Among 
these parameters, J, = 160560 kg - m’, and Di=0 = 
1.7827 x 10° N - m - (rad - s7)! are already known (Chen 
& Wang 2022); and D2, Kz can be calculated from the 
following equations (Gawronski 2008): 


Dz = 26, W2Ja (8) 
Ky = w7*Jn (9) 


where ¢,;=0.01 is the damping ratio of the main structural 
vibration modal; and the angular velocity w2 of the main 
structural vibration modal is 


w = Inf, = 2 x m x 13.004 = 81.7065 rad : s! (10) 


(6) 


where f> = 13.004 Hz is the natural frequency of the structural 
vibration modal (1.e., the second-order vibration modal of the 
antenna). 
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By applying Laplace transform on both sides of 
Equations (6) and (7), the input-output model of LCT’s 
antenna can be obtained as below: 


ro 
Jas? + Dis 
(11) 
ho) = ——) 
Jas? ae Das + Kə 


By integrating the model of LCT’s antenna with two major 
vibration modals to the pre-established LCT’s pointing control 
system model (Chen & Wang 2022), we can obtain the 
extended open-loop model of the pointing control system as 
shown in Figure 1. 

Note that Tz = Ts if the wind disturbance is not introduced to 
the system, where T, is the output torque of the DC motor and 
reducer. Assume that xı =aRM, %2=Wpcm, X3 = ĝÎpcCM, 
x4 = 1, X5 = Way, Xo = Oo, X7 = War, Where Iapy is the armature 
current of the DC motor; wpcm is the rotating speed of the DC 
motor; wa; and wy are the rotating speed of the rigid body of 
the antenna, and the rotating speed of the main structural 
vibration modal of the antenna respectively; Opcm is the 
rotation angle of the DC motor. Then the state-space model of 
LCT’s pointing control system with extended LCT antenna 
model can be written as follows: 


_ Uxrm/Ucom Rpcm Cocm 
Uc x X2 
Lpcm Lpcm Lpcm 
1 
x2 = Kpcm™ — o1X2 — En) 
ml N 
x3 = x2 
X4 = x5 (12) 
D 1 
ïs = -xs + —Ty 
Ja Ja 
X6 = X7 
K D> 1 
X% = X6 X | Tz 
Ja Ja Ja 
180 
y = — (x4 + x6) (13) 
T 
with 
K 
i= x — Kepc(x4 + x6) (14) 


where x; (i = 1...7) are the system states, y is the system output 
(i.e., the position of the antenna in azimuth, or the azimuth 
angle, in degrees), u, is the control input (i.e., the command 
voltage with the maximum value of Ucom); and the parameters 
of Equations (12)—(14) are listed in Table 1. Note that T; in 
Equations (12) and (14) is the load torque of the reducer, and 
T,/N is the load torque of the DC motor. 

It should be noted that the value of Kpcm (i.e., the torque 
coefficient of the DC motor) in Table 1 is twice the value of the 


Chen & Wang 


real single DC motor, because there are two DC motors in the 
azimuth channel of the real-world LCT’s pointing control 
system for the purpose of eliminating the backlash of the gears 
in the speed reduction system and they are equivalent to one 
DC motor in this model. 

From Equations (12)-(14), the state-space model of LCT’s 
pointing control system with extended LCT antenna model can 
be formulated as below: 


f = AsX; + Bye 


y= CyXs + Dye uo) 


where x, = [x x2 X3 x4 X5 Xo X7]" is the state vector; A,, Bs, Co 
D, are the system’s parameter matrices that can be formulated 
as follows: 


— Roem _ Coc 0 0 0 0 0 
Lycm Loom 
Kpcm ai Kroc Kroc 0 Kroc 0 
Jocm Jocm  JocmN? JpcmN JocemN 
0 1 0 0 0 0 0 
A= 0 0 0 0 1 0 0 
0 0 Kroc Kroc Dı Kroc 0 
JaN Ja Ja Ja 
0 0 0 0 1 0 0 
0 0 Kroc Kroc 0 Kroc +K Da 
JaN Ja J Ja 
(16) 
Uarm/ Ucom i 
B = | —— 000000 (17) 
Lpcm 


C=[0 0 0 180/r 0 180/r 0],D=0 (18) 


3. Design of the DOB-Based Position Controller 


There exist two types of disturbances in LCT’s pointing 
control system: one is the external disturbance, i.e., the wind 
load, and the other is the internal disturbance, i.e., the vibration 
of the flexible structure of the antenna. The two types of 
disturbances are coupled with each other when they act 
simultaneously. Whereas the DOB is known for its good 
performance on estimating the external and internal disturbance 
accurately, it is expected that the negative impacts of the 
external and internal disturbance can be mostly rejected if a 
well-designed DOB-based position controller (i.e., a position 
controller composed of a DOB-based feedforward controller 
and a feedback controller) is employed in LCT’s pointing 
control system. Hence, in this section, we focus on the design 
of the DOB-based controller of LCT’s pointing control system. 
First, we design a time-domain DOB to estimate the coupled 
disturbance generated by both the external disturbance (i.e., the 
wind load) and the internal disturbance (i.e., the structural 
vibration modals). Second, a suitable linear-quadratic-regulator 
(LQR) - based feedforward controller is designed to compen- 
sate the disturbance estimated by the DOB. Finally, we design 
a composite DOB-based position controller by integrating the 
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Figure 1. Diagram of the extended open-loop model of LCT’s pointing control system with structural vibration modals. 
Table 1 
Parameters of LCT’s Pointing Control System 
Symbol Definition Value 
Ja the moment of inertia of the LCT antenna 160560 kg - m? 
D; the damping of the rigid body of the LCT antenna 1.7827 x 106 N - m- (rad - s~!) ! 
D2 the damping of the main structural vibration modal 2.6237 x 10° N - m - (rad - s7}! 
Kə the stiffness of the main structural vibration modal 1.07 x 10° N - m - rad™’ 
Kerpc the stiffness of the reducer 3.498 x 107N-m-rad ! 
Lpcm the inductance of the DC motor 0.00159 H 
Rpcm the resistance of the DC motor 0.30 ohm 
Kpcm the torque coefficient of the DC motor 0.542 N- m- A7! 
Jpcm the moment of inertia of the DC motor 2.034 x 10 7N-m-s 
Cpcm the back electromotive force (EMF) coefficient of the DC motor 0.271 V - (rad - s7!) 
o the damping coefficient of the DC motor 4.8694 x 107° N - m - (rad - s7)! 
N the reduction ratio 15 050 
UARM the armature voltage of the DC motor 115 V 
Ucom 


the maximum value of the command voltage 5V 


feedforward controller and the feedback controller (Ha 
controller or proportional controller), along with a controller 
switching policy for keeping good transient and steady-state 
performance while rejecting external and internal disturbance. 


3.1. The Time Domain DOB (TD-DOB) 


There exist two types of the DOB, which are the frequency 
domain DOB (FD-DOB) and the time domain DOB (TD-DOB) 
(Li et al. 2016). By taking the antenna as the controlled object 
(i.e., the object to be controlled by the control system) for more 
understandable explanations, the key idea of FD-DOB can be 
explained as follows: If we know the nominal model of the 
antenna, then we can estimate the nominal input torque signal 
from the inverse of the transfer function of the antenna’s rigid 
model and the measured output azimuth angle of the antenna. 
In the absence of external disturbance (i.e. the wind 
disturbance) and the internal disturbance (1.e., the vibration of 
the flexible structure of the antenna), this estimated input torque 
signal should be consistent with the output torque signal of the 
motor and reducer. However, in practice, the input torque of the 
antenna is affected by the torque of the external disturbance 


(such as the wind disturbance), and the antenna’s flexible 
modals lead to the fluctuation of the output azimuth angle, 
which makes the actual input torque not equal to the nominal 
input torque. Then, the input torque disturbances caused by 
internal disturbance (i.e., the vibration of the flexible structure 
of the antenna) and external disturbance (i.e., the wind 
disturbance) can be finally obtained by differentiating the 
actual output torque from the calculated nominal input torque. 
By estimating the input torque disturbances and compensating 
it accordingly by regulating the motor, the negative effects of 
the input torque disturbances on the pointing control system 
can be counteracted to some extent. 

TD-DOB is another formulation of the DOB. Similar to FD- 
DOB, TD-DOB is motivated by the desire of using the model 
information and the input and output signals of the model to 
extract the disturbance signals from the input signals. There- 
fore, TD-DOB also requires a nominal model of the antenna, 
the measured output azimuth angle, and the output torque 
signal from the motor and reducer. However, unlike FD-DOB 
which directly utilizes the inverse of the transfer function of the 
antenna’s rigid model to estimate the disturbance signal in the 
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frequency domain, TD-DOB is constructed in the time domain 
by augmenting the antenna model in the form of a state-space 
by introducing a disturbance term, and thus constructing an 
observer to estimate the augmented state including the 
disturbances. 

Considering that the TD-DOB is convenient to implement on 
embedded hardware platforms, we adopt the TD-DOB for 
disturbance estimation. At first, to estimate the composite 
disturbance generated by both the external disturbance (i.e., the 
wind load) and the internal disturbance (i.e., the structural 
vibration modals), we take the rigid body model of LCT’s 
antenna as a nominal one, which can be obtained from 
Equation (12) as the following state-space equations: 


i = AnXa + B, (T; + Ta) 


19 
0a = CiXa “ 


where x, = [war 0] is the vector of the states of the rigid 
body model of the antenna; T4 is the torque generated by the 
wind load; and A,, B,, C, are the matrices of the system 
parameters defined as below: 


-Pi o 2 
An = Ja , B, = JA |> C = [o 1] (20) 
1 0 0 


To facilitate the implementation of the DOB on embedded 
hardware platforms, inspired by Yang et al. (2017), the TD- 
DOB is designed as follows 

Z= LB, (Zz t Lxa) 
Ù =z + Lya 


L AnXa + BrT 
( ) (21) 


where 7, is the estimated torque of the internal and external 
disturbance; L > 0 is the gain matrix of the TD-DOB, which 
needs to be further designed; z is the auxiliary variable of the 
DOB; and T, is the nominal load torque of LCT’s antenna 
without any disturbances. It should be noted that 7, can hardly 
be measured directly and hence it must be calculated according 
to the position of the antenna 04 = 0; + 02 and the speed of the 
motor Wpcm as follows: 


RO = Keoe( Z f evca(riar - a0) (22) 


Moreover, a first-order low-pass filter is designed to provide 
the proper estimated disturbance signal for designing a 
realizable feedforward controller, which is given by 


P(s) = 


23 
As +1 a 


where A € R* is the parameter of the first-order low-pass filter. 
Therefore, the ultimate output signal To of the TD-DOB is 


p= P= (24) 
As +1 
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Figure 2. Diagram of the TD-DOB for LCT’s pointing control system. 


By combining Equations (19)-(24), the overall structure of 
the TD-DOB is designed as Figure 2 shows. The designed TD- 
DOB has the capability of estimating the composite external 
and internal disturbances since the output of the whole antenna 
model, rather than the output of the rigid body part of the 
antenna’s model, is input to the TD-DOB. 


3.2. LOR-Based Feedforward Controller for the DOB 
(LQR-DOB) 


The output of the TD-DOB (..e., Tip in Equation (24)) is an 
estimation of the torque generated by the composite external 
and internal disturbances, which will be used to design a 
feedforward controller to compensate the real disturbance. A 
straight-forward idea is that if we can design a feedforward 
controller such that the output torque of the actuator (which 
includes the driver, the DC motors, and the reducer) can track 
Tio accurately, then the torque generated by the disturbance will 
be compensated quite well. 

According to LCT’s pointing control system captured by 
Equation (15), we construct the state-space model that 
represents the relationship between the feedforward control 
signal (which is denoted by ug) and the output torque of the 
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actuator (which is denoted by T,) as follows: 


Xq = AgXq + Betta 
= (25) 
Ta = CyXa 


h — [ya d a T is th p 
[where x4 = IRM “pcm ODcM 1s the state vector, 1n 


which fkm wfcm and O44 are the armature current, the 
rotating speed, and the rotating angle of the DC motor 
respectively; and the matrices A,, B,, C, are formulated as 
follows, 


_Roem _ Coe o 
Lpcm = Lpcm 


Ag=| Kom __% oF (26) 
Jpcm JpcM 
0 1 0 
x= | Uarm/ Ucom 0 o], (27) 
Lpcm 
Kroc 
&=fo o Kew | 28) 


In order to make 7; accurately track the output of the DOB, 
i.e., Zip, the estimation of the torque generated by the composite 
external and internal disturbances, a linear quadratic optimal 
tracking control problem is formulated as follows: 


min J = =f [eT (t) Oe(t) + uj (t)Rua(t)] (29) 
ža(t) = Agxa (t) + Byua(t), Xa(to) = x} 

s.t.) Ty(t) = Cyxa(t) (30) 
e(t) = Tyo(t) — Ty(t) 


where e(f) is the tracking error; Q € R is the weight of the 
tracking error; R € R is the weight of the feedforward control 
signal; and x? is the initial state of the system. 

According to the optimal control theory (Lewis et al. 2012), 
a linear quadratic regulator (LQR) can be designed to generate 
the approximate optimal control signal as follows 


ug(t) = —R'B,7Pxq(t) + R7'B,78 (31) 


where P is a positive definite matrix satisfying the following 
Algebraic Riccati Equation (ARE) 


PA, + ATP — PB,R—B,'P + C,™@C, = 0 (32) 
with 
ê = [PB,R-'B,? — A,T'C,7QTio (33) 


Define G = A, — B,R7'B, TP then the LQR-based feedfor- 
ward controller for the DOB (LQR-DOB) can be finally 
designed as Figure 3 shows, where the estimated disturbance 
signal To) is the input of the feedforward controller and u(t) 
is the output of the feedforward controller. Note that the 
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amplitude of the control signal u4(t) can be limited by adjusting 
the parameters Q and R. 


3.3. DOB-Based Composite Position Controller with the 
Controller Switching Policy (DOB-CPC with CSP) 


To reject the disturbance from LCT’s pointing control 
system, the TD-DOB, LQR-based feedforward controller, and 
the feedback controller are combined as the DOB-based 
composite position controller (DOB-CPC) to replace the 
existing position controller (which is a proportional controller). 
Therefore, the control signal u. generated by DOB-CPC 
consists of two parts: the feedback control signal u, and the 
feedforward control signal ug. 

The first part is the feedback control signal u, from the 
feedback controller (i.e., the existing proportional feedback 
controller or the H,, feedback controller proposed in our 
previous work (Chen & Wang 2022)), which calculates the 
control signal based on the position error of LCT’s antenna. As 
LCT’s pointing control system is captured by Equation (15), 
the Ha feedback controller can be obtained in many ways. 
Here the Ha loop-shaping technique (Glover & McFarlane 
1989; Sefton & Glover 1990; Feyel 2013) is used to design the 
H feedback controller based on the model of LCT’s pointing 
control system because it is easy to implement by designing the 
three loop-shaping weighting filters W1, W2, W3 with the help 
of the MATLAB Robust Control Toolbox.* As Figure 4 shows, 
the error signal e, the control signal u, and the output signal 
y are input to the loop-shaping weighting filters W,, Wo, 
W3 (which can be set as the constant, the high-pass filter, 
or the low-pass filter) respectively to generate the measures of 
robustness z1, Z2, Z3 (which are parts of the infinity norm of the 
closed-loop transfer function of the augmented model of LCT’s 
pointing control system). Since the design of the Ha feedback 
controller aims to minimize the infinity norm of the closed loop 
transfer function of the pointing control system, the error signal 
e, the control signal uc, and the output signal y can be adjusted 
by tuning Wi, W2, W3 respectively. 

The second part is the feedforward control signal u4 from the 
LQR-based feedforward controller, which is obtained based on 
the estimated disturbance signal i. by using TD-DOB. Hence, 
the DOB-CPC not only makes full use of the characteristics of 
the feedback controller to quickly eliminate the position error, 
but also compensates the disturbance to further reduce the 
position fluctuation it generates. 

Furthermore, as the step response plots shown in Figure 5, it 
should be noted that the Ha feedback controller and the LQR- 
based feedforward controller will bring negative impact on the 
rapidity of the transient process although they can greatly 
enhance the ability of disturbance rejection, especially for the 
circumstances under high-speed wind disturbances. 


4 https: //www.mathworks.com/help/releases /R2022a/robust/ref/Iti. 
hinfsyn.html 
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Figure 3. Diagram of the LQR-based feedforward controller for the DOB. 


In specific, in the absence of wind disturbances (See 
Figure 5(a)), the transient process to the pointing control 
system with the H,, feedback controller is longer compared to 
the proportional (P) controller while the LQR-based feedfor- 
ward controller (LQR-DOB) does not affect the pointing 
control system’s transient performance. However, this situation 
changes significantly when facing large wind disturbances (See 
Figure 5(b)). In addition to the negative impact of the Ha 
feedback controller on the pointing control system’s transient 
performance, the LQR-based feedforward controller (LQR- 
DOB) often generates overcompensation in the transient 
process, which leads to a significant delay in response. Thus, 
considering that the transient process under the existing 
proportional controller is the fastest (since the control signal 
is maintained at the maximum level), we propose a controller 
switching policy (CSP) for DOB-CPC. The key idea of the 
CSP is to activate the H,, feedback controller and the LQR- 
based feedforward controller only when the position of the 
antenna is within a certain range around the desired position 
(e.g., the orange area shown in Figure 6); otherwise, the 
existing proportional feedback controller works. 

For the purpose of utilizing different controllers for the 
transient and steady-state stages respectively, then the CSP has 
the following form: 

f = up, |e,| > Ah 


34 
Uc = kpup + kjua, |e,| < Ah a 


where e, is the position error of the antenna; Ah € [0, r] is the 
threshold of position error for activating H,, feedback control 
(which can also be replaced by the existing proportional 
feedback control, namely, keep using the existing proportional 
feedback control) and feedforward control; u, is the feedback 
control signal; u, is the control signal obtained by the existing 
proportional feedback controller; ug is the control signal 
obtained by the LQR-based feedforward controller using the 
DOB-estimated disturbance signal; kp € R is the weight of the 
feedback control; and ky € R is the weight of the feedforward 
control. Here the equation that satisfies |e,| > Ah is for the 
transient process, and the equation that satisfies |e,| < Ah is 


Pointing Control 
System Model 


Controller 


a 


Figure 4. The diagram of closed-loop LCT’s pointing control system with Hao 
feedback controller and loop-shaping weighting filters W,, W2, W3: r is the 
reference input, e is the error signal, uo is the control signal calculated by the 
H feedback controller, y is the system output signal, z1, z2, z3 are the measures 
of robustness. 


for the steady-state process. By employing the CSP, the 
structure of DOB-CPC is depicted as the blue block in 
Figure 7. 


4. Simulation Experiments 
4.1. Experiment Design 


In this section, the experiment is divided into two parts: 

(1) Mathematical Simulation Experiment. MATLAB/ 
Simulink is used to perform the simulation experiments for 
testing the performance of the proposed composite DOB-based 
position controller based on the mathematical model of LCT’s 
pointing control system (whose antenna model includes a rigid 
body and the main structural modal). In this experiment, the 
performance of four types of position controllers (which are 
listed in Table 2) under the wind disturbance will be tested and 
compared. 

(2) Collaborative Simulation Experiment: Due to the 
unavailability of LCT during its relocation process, we 
developed a collaborative simulation testbed (Yao et al. 
2023) to verify the performance of the proposed DOB-CPC 
(H,.+LQR-DOB) by introducing a pre-constructed multibody 
dynamic model of the antenna. 
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Figure 5. Step responses of the pointing control system using different controllers for two situations: (a) without wind disturbances; and (b) with 20 m- s~' wind 
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Figure 6. An example of the position of LCT’s antenna. The range for 
activating the H,, feedback controller and the LQR-based feedforward 
controller is marked as the orange area. 


As described in our previous research article (Yao et al. 
2023), the Adams model is constructed based on the Ansys 
model of LCT. To be specified, as Figure 8 shows, we first 
constructed the mechanical structure of LCT using SpaceClaim 
(a 3D modeling tool within Ansys). The bodies of the 
mechanical structure have two categories, i.e., rigid bodies 
and flexible bodies. On one hand, rigid bodies can be well 
supported by Adams so that we can import them to Adams 
directly. On the other hand, flexible bodies are difficult to 
import to Adams directly since Adams is not good at building 
complex structural bodies. Therefore, we utilize Ansys to 
transform LCT’s mechanical structure model to a finite 
elements model and proceed the modal analysis of LCT based 
on its finite elements model to obtain the modal neutral file of 
the flexible bodies. Finally, with the help of Adams, the rigid 
bodies (imported directly) and the flexible bodies (imported by 
the modal neural file) are integrated as a whole by applying the 
fixed joints to the contact surfaces between the rigid bodies and 
the flexible bodies. 

With the help of the established Adams model, the 
simulation testbed is able to reflect the complex multimodal 
characteristics of the antenna in Simulink and its impacts on the 


dynamics of the entire pointing control system. (Without 
Adams, the complex multimodal characteristics of the antenna 
can only be analyzed statically in Ansys, which was not able to 
analyze the influence of these modals on the dynamics of the 
entire pointing control system). 

Based on the existing model of the pointing control system, 
we replace the state-space formulation of the antenna by its 
multibody dynamics model, which was constructed using 
Adams (a multibody dynamics simulation software) by our 
team (Yao et al. 2023) to reduce the gap between the simulation 
model and the antenna in the real world. In addition, the 
multibody dynamics model of LCT’s antenna has the same 
properties (i.e., material, dimensions, structures, density, and 
stiffness) as the antenna in the real world according to the on- 
site measurement and the Ansys model provided by our 
Chilean colleagues (See Acknowledgments) such that the 
accuracy of the multibody dynamics model of the antenna can 
be guaranteed. Then, with the collaborative simulation 
combining MATLAB /Simulink and Adams, we can further 
verify the proposed DOB-CPC (H,.+LQR-DOB) on the 
collaborative simulation testbed (i.e., an integration of the 
mathematical model and the multibody dynamics model of the 
antenna, as shown in Figure 9). In this experiment, the 
proposed DOB-CPC (H,.+LQR-DOB) will be tested on the 
mathematical model and the collaborative simulation testbed 
under different wind disturbances to observe whether it still 
works well in the situation close to the real world. 

The wind disturbance consists of the wind gust (captured by 
the fluctuation of the wind speed, i.e., Av,,) and the steady 
wind (captured by the average speed of the wind, i.e., v,,), in 
which the wind gust is the main factor leading to the rise of 
position error. Therefore, it is important to construct a good 
model for wind gust. To better describe wind gust, the 
Davenport filter is introduced, which is designed by following 
the instructions in Section 3 of Chen & Wang (2022). The 
power of the white noise (i.e., the source of the Davenport 
filter) and the average wind speed are adjustable for generating 
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Figure 7. Diagram of LCT’s pointing control system model when the DOB-based composite position controller (DOB-CPC) with CSP is implemented. “P” refers to 
the existing proportional feedback controller. “Hao” refers to the Ha feedback controller. 
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Figure 8. Diagram of the construction of the Adams model of LCT. 


Table 2 
Descriptions of Different Position Controllers to be Tested 


Controller Type 


P 
Ho 
DOB-CPC (P+LQR-DOB) 


Description 


The H» feedback controller 


The existing proportional feedback controller 


The DOB-based composite controller (DOB-CPC) shown in Figure 7 activates the existing proportional feedback controller and the 


LQR-based feedforward controller for the DOB (LQR-DOB) on the steady-state stage 


DOB-CPC (H,,+LQR-DOB) 


The DOB-based composite controller (DOB-CPC) shown in Figure 7 that activating the Ha feedback controller and the LQR-based 


feedforward controller for the DOB (LQR-DOB) on the steady-state stage 


wind disturbances of different average speeds and fluctuation 
amplitudes. According to the detail of the climate condition of 
the new site at Chajnantor Plateau, Chile, on one hand, the 
wind speed changes slowly during a day; on the other hand, the 
maximum amplitude of the wind speed fluctuation is 5 ms ' 
(i.e., the wind speed fluctuates from 6 to 11 ms! during a day 
in winter) (ALMA 2010). Additionally, the average wind speed 
of the new site is 10 m s |, and the maximum wind speed of 
the new site is 30 ms ' (ALMA 2010). Based on the above 
information, we have tested and compared the position control 
performance of four types of controllers under different average 
wind speeds and different wind gust fluctuations, where the 
average wind speed is set to three levels of 10, 15, and 
20 ms_', and the amplitude of wind gust fluctuation is set to 
three levels of 2, 6, and 10 ms !. For convenience, the 
amplitude of wind gust fluctuations is measured by the 


10 


following equation (see Figure 10): 


Aay,, = |Avycmax) = AVw(min)| (35) 


where Aav, is the amplitude of wind gust fluctuations, Avw(max) 
is the maximum value of the wind gust fluctuations, and 
AYVw(min) is the minimum value of the wind gust fluctuations. 
Then, the wind speed is converted to the equivalent wind 
disturbance torque (T4) applied to the antenna’s azimuth axis 
by multiplying a conversion factor (which is calculated based 
on the wind area of the antenna, the static air density, and the 
wind direction) according to the method described in (Chen & 
Wang 2022). 

For the fairness of comparison, all controllers will be well 
tuned and then be tested in the same situation (i.e., the same 
target position (i.e., the azimuth angle) of the antenna, the same 
pointing control system model, and the same wind disturbance). 
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Figure 9. Collaborative simulation testbed for LCT’s pointing control system (which has the Adams model of the antenna) when a DOB-based composite position 
controller (H,,+LQR-DOB) is implemented. o4 = 1.7827 x 106 N -m - (rad - s~!)~! is the damping of LCT antenna. 


We observe the output of LCT’s pointing control system model ( 
i.e., the position of the antenna) with different controllers (i.e., 
the existing proportional feedback controller, the Ha feedback 
controller, the DOB-CPC with existing proportional feedback 
controller, and DOB-CPC with H,, feedback controller) under 
the same wind disturbance, when the target position (i.e., the 
azimuth angle) is set as 10 deg. 

The controllers are pre-configured as follows: First, the 
formulation and parameters of the existing proportional 
feedback controller are kept as given by the original technical 
report (Leighton 1977). Second, the H,, feedback controller is 
redesigned based on the state-space model of the pointing 
control system with vibration modals (see Equations (12))— 
(14)) by using the following loop-shaping weighting filters Wj, 
W2, W3: (1) W, is set as a low-pass filter such that the low- 
frequency gain (LFG) of W; is 5000 (i.e., 74 dB), the frequency 
corresponding to 0dB magnitude of W; is 5 rad - s_', and the 
high-frequency gain (HFG) of W; is 0.1 (i.e., —20 dB); (2) W2 
is set as a constant, which is equivalent to 0.000 001; (3) W3 is 
set as a high-pass filter such that the LFG of W3 is 0.1 (ie., 
—20 dB), the frequency corresponding to 0dB magnitude is 
6 rad - s7}, and the HFG of W; is 5000 (i.e., 74 dB). The Bode 
diagram of W, W3 is illustrated in Figure 11. Third, the 
parameters of the proportional feedback controller and the Ha 
feedback controller used in the DOB-CPC are the same as 
stated in the first two points, and the other undeclared 
parameters of DOB-CPC are tuned one by one for each case. 


4.2. Simulation Results and Discussions 


A. Mathematical Simulation Experiment 

To carry out the mathematical simulation experiments, the 
parameters of the controllers are first well tuned, including: (1) 
the TD-DOB’s parameters L, A, which are the gain matrix of 
the TD-DOB and the parameter of the first-order low-pass filter 
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Figure 10. Diagram of the wind gust fluctuations. 


P(s) of the TD-DOB respectively (which can be found in 
Equations (21) and (23)); (2) the LQR-DOB’s parameters Q, R, 
which are the weights of the tracking error and the feedforward 
control signal of the LQR-DOB respectively (which can be 
found in Equation (29)); (3) the CSP’s parameters kp, ky, Ah, 
which are the weight of the feedback control signal u, 
calculated by the feedback controller (proportional or H), 
the weight of the feedforward control signal ug calculated by 
the LQR-DOB, and the threshold of position error for 
activating Hə feedback control (or the existing proportional 
feedback control) and feedforward control respectively (which 
can be found in Equation (34)). Part of the parameters of the 
controllers are configured as Table 3 shows, and other unlisted 
parameters are set as follows: Q=15, k,=1, ky=1 and 
Ah=0.1 are adopted for all cases of simulation experiments. 

Figures 12—14 present the position error of LCT’s antenna 
using different controllers under different average wind speeds 
and the fluctuation amplitudes of 2, 6, and 10 m s 
respectively. To further quantify and measure the disturbance 
rejection performance of the controllers, the root mean square 
(rms) of the position error is measured during the interval from 
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Figure 11. Bode diagram of loop-shaping filters W), W3. 


the 10th s to the 20th s. Table 4 concludes the RMSs of the 
position error of LCT using different controllers under different 
average wind speeds and the fluctuation amplitudes of 2, 6, and 
10 ms™!. 

From Figures 12-14, we can see that the transient 
performances of the pointing control system under the four 
controllers are almost the same because of the effectiveness of 
the CSP, i.e., the same proportional control is adopted before 
entering the steady state. Then, when the position error is 
within the range for activating Hə control or feedforward 
control, the switching of controllers begins to show their 
strength on wind disturbance rejection. From Figures 12—14 
and Table 4, we can see that the employment of the DOB-CPC 
controller (especially the H,.+LQR-DOB controller) can 
reduce the rms of the position error significantly. In addition, 
Figures 12—14 and Table 4 depict a clear trend that the position 
error increases with the rising of the wind average speed and 
the wind fluctuation amplitudes. However, even if under a high 
wind speed of 20 ms’ and a large wind fluctuation of 
10 ms’ |, the rms of position error can be reduced by more than 
90%. For a lower wind speed of 10 ms_' and a smaller wind 
fluctuation of 2 ms ', the rms of position error can be reduced 
by more than 92%. 

On one hand, by comparing the position error under the 
existing proportional controller, the H,, feedback controller, 
and the DOB-CPC (P+LQR-DOB), the position error of the 
LCT’s pointing control system controlled by the Ha feedback 
controller is first found much smaller than the one controlled by 
the existing proportional feedback controller. Moreover, once 
the DOB-CPC (P+LQR-DOB) is employed, the position error 
of the antenna is dramatically reduced compared to the one 
under the Ha controller or the existing proportional feedback 
controller. Obviously, replacing the existing proportional 
controller by the H, controller, or combining the existing 
proportional controller with the LQR-DOB feedforward 


12 


Chen & Wang 


Table 3 
Part of the Parameters of the Controllers for the Wind Disturbance of Different 
Average Wind Speeds and Fluctuation Amplitudes 


Case Description iL R AÀ 

vw = 10m s~!, Aan, = 6ms™! [10074 100J,] 0.3 0.001 
vp = 10ms7!, Aan, = 10 m s7! [10074 100J,] 0.3 0.001 
w = 15 m s™!, Aay =2 ms! [10074 100J,] 0.3 0.001 
ww = 15 m s7!, Aan, = 6 ms! [10074 100J,] 0.3 0.001 
ww = 15 m s™!, Aay, = 10m s7! [10074 100J, 0.3 0.001 
vy = 20ms"!, Aan, = 2 ms! [10074 100J, 0.3 0.001 
vw = 20ms—!, Aan, = 6m s7! [10074 100J, 0.3 0.003 
ww = 20 m s7!, Aav, = 10 m s7! [10044 100J4 0.3 0.003 


controller (i.e., the DOB-CPC controller), is helpful in reducing 
the position error of the antenna under the wind disturbance. In 
addition, the introduction of LQR-DOB feedforward controller 
contributes much more to eliminating the negative effect of 
wind disturbance on position control, compared with the 
contribution of the H,, controller. This phenomenon can be 
explained by the principles of the Hə controller and the LQR- 
DOB as follows: (1) The H, controller is capable of 
compressing the range of the position error to reduce the 
position error globally since it can improve the robustness of 
the system. (2) The LQR-DOB (which is integrated with the 
feedback controller as the DOB-CPC) provides an option of 
eliminating the negative effect of the external and internal 
disturbances precisely by estimating the disturbance signal 
using the DOB. Therefore, it achieves better performance 
compared to the Hæ controller. 

On the other hand, by comparing the results of using the DOB- 
CPC (P+LQR-DOB) and the DOB-CPC (H,,+LQR-DOB), it 
can be found that the position error of the antenna is further 
reduced when the feedback controller of the DOB-CPC is 
changed from the existing proportional feedback controller to the 
Ha feedback controller. A possible explanation for this 
phenomenon is that the combination of the Ha controller and 
the LQR-DOB feedforward controller combines the strengths of 
both: rough disturbance rejection (compression of the range of 
position errors) with the Hə» feedback controller, and more 
refined disturbance rejection with the LQR-DOB feedforward 
controller. 

The pointing accuracy of LCT is 3 arcsecond (3”) (rms), 
which requires that the rms of the position error of its pointing 
control must be less than 3”. However, the RMSs of the 
position errors in all cases using the existing proportional 
feedback controller exceed 3”, as shown in Table 4. After 
applying the proposed DOB-CPC (H,,+LQR-DOB), the 
RMSs of the position errors in all cases are greatly reduced 
by more than 90% accompanied by most cases achieving less 
than 3”. Nevertheless, the proposed DOB-CPC may have some 
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Figure 12. Position error of LCT’s antenna using different controllers under different average wind speeds and the fluctuation amplitude of 2 m s™'. (a) 
vw = 10ms"!, Aaw, = 2ms! (b) vw = 15 m s7}, Aav, =2ms! (c) vw = 20ms"!, Aav, = 2 m s7! 
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Figure 13. Position error of LCT’s antenna using different controllers under different average wind speeds and the fluctuation amplitude of 6 m s™'. (a) 
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Figure 14. Position error of LCT’s antenna using different controllers under different average wind speeds and the fluctuation amplitude of 10 ms ~. 
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Table 4 
RMSs of the Position Errors of LCT’s Antenna Using Different Controllers 

DOB-CPC DOB-CPC % of rms 
Case Description P He (P+LQR-DOB) (Hao + LQR-DOB) Decrease 
vw = 10 ms™!, Aav, = 2 ms“! 4626 3/173 0/461 0/366 92.09% 
vw = 10 ms~!, Aav, = 6 ms! 13”860 9534 1/377 1001 92.78% 
vy = 10 ms~!, Aav, = 10 ms“! 237110 15”900 27348 1680 92.73% 
vw = 15 ms}, Aay, = 2 ms“! 87946 6/108 0/872 0/648 92.76% 
vp = 15ms~!, Aav, = 6 ms! 267770 187320 21800 1/987 92.58% 
vy = 15 ms, Aav, = 10 ms“! 44800 30/670 6/466 4479 90.00% 
vy = 20ms~!, Aay, = 2 ms! 16”850 117590 1/673 1/214 92.80% 
Vw = 20 ms~!, Aav, = 6 ms“! 507540 347800 51856 4/230 91.63% 
Yw = 20 ms~!, Aav, = 10 ms“! 84210 58/000 10/920 7!'856 90.67% 


Note. Under Different Wind Disturbances of Different Average Speeds and Fluctuation Amplitudes. 


limitations on its performance for some extreme cases, which 
are (1) the wind disturbances with an average speed of 
15 ms! and the fluctuation of 10 ms! and (2) the wind 
disturbances with the average speed of 20 ms’ and the 
fluctuation above 6 ms '. The above results provide some 
instructions of applying the DOB-CPC properly. If a high 
pointing accuracy (e.g., 3”) is required during LCT’s operation, 
employing the DOB-CPC can achieve this aim when the wind 
disturbance is not very large (e.g., a relatively low average 
speed of 10 ms_', or a relatively small fluctuation of 2 m s7». 
If the wind disturbance is relatively large (i.e., an average speed 
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of 20 m s |, or a relative medium fluctuation of 6 m s7!) but 
the requirement on the pointing accuracy is not very high (e.g., 
above 3” but below 10”), the DOB-CPC still works well. 

In general, by carrying out the above tests, we find that the 
proposed DOB-CPC can help LCT’s pointing control system to 
deal with different types of wind disturbances and achieve 
significant performance improvement compared to the single 
feedback controller (i.e., the proportional controller or the Ha 
controller). Thus, the accurate estimation of disturbance signals 
by TD-DOB and the accurate generation of the feedforward 
control signal by LQR together eliminate the negative effect of 
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Figure 15. Position error of LCT antenna using the proposed DOB-CPC (H,,+LQR-DOB) on both the mathematical model and the collaborative simulation testbed 
under different wind disturbances (a) v = 10 m s™!, Aaw, = 2 m s™! (b) v = 15 m s™!, Aav, = 6 m s™! (c) vy = 20 m s™!, Aaw, = 10 m s™!. 


the external and internal disturbances caused by wind load and 
structural flexibility of the antenna. 

B. Collaborative Simulation Experiment 

By integrating LCT’s pointing control system model and the 
multibody dynamics model of LCT’s antenna, we can conduct 
the collaborative simulation experiment to further verify the 
performance of the proposed DOB-CPC (H,.+-LQR-DOB) in the 
situation closer to reality. This experiment is conducted under 
three typical wind disturbances, which are (a) v,= 10ms_}, 
Aq =2 ms! (b) w=15ms™}, Ag Homes (c) w= 
20ms™!, Aay, = 10 m s™!. Considering that the dynamics of 
the antenna in the state-space formulation is different from the 
dynamics of the antenna in the Adams model (which contains 
more vibration modals and nonlinear characteristics that are not 
captured by the state-space formulation), the former well-tuned 
parameters of the controllers need to be adjusted for adapting the 
Adams model in the collaborative simulation testbed: (1) change 
AÀ from 0.001 to 0.006 for the case of v, = 10ms—!, 
Aaw, = 2ms™!; (2) change » from v, = 15 m s™!, Aa, = 
6m s™!; (3) change from 0.003 to 0.02 for the case of 
vy = 20 m s™!, Aay, =10 ms}. 

Figure 15 shows the position error of LCT’s antenna under 
the DOB-CPC (H,.+LQR-DOB) for both the mathematical 
model (i.e., the state space model) and the collaborative 
simulation testbed (i.e., the state space model integrated with 
the Adams model of antenna) under different wind distur- 
bances, and the RMSs of the position errors are also calculated 
and shown in Table 5. Figure 15 and Table 5 indicate that the 
performance of the proposed DOB-CPC (H,.+LQR-DOB) on 
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the collaborative simulation testbed is almost the same as the 
one on the mathematical model under the most often-seen wind 
disturbance (i.e., vy = 10 m s™!, Aay, = 2ms~'). However, 
for stronger wind disturbances, the proposed DOB-CPC 
(H,.+LQR-DOB) performs relatively worse on the collabora- 
tive simulation testbed than on the mathematical model: the 
position error increases from 1”9860 to 4”0290 for the case of 
vw = 15 m s™!, Aay, = 6 m s™!, and from 78560 to 18/3300 
for the case of v, = 20 m s™!, Aay, = 10 m s™!. Despite of the 
degradation of performance in these two cases, the position 
error of 4”0290 is still acceptable under the average wind speed 
of 15 ms™' and the fluctuation of 6 ms’ '. Actually, at LCT’s 
new site on Chajnantor Plateau in Chile, the first two cases of 
wind load are much more common than the third case. 
Therefore, the result of this experiment verifies the feasibility 
of applying the proposed DOB-CPC (H,,+LQR-DOB) on 
LCT’s pointing control system after it is relocated and 
refurbished in Chile. 


5. Conclusions 


In this paper, a disturbance observer-based composite 
position controller (DOB-CPC) for LCT’s pointing control 
system is proposed to realize a better performance of 
disturbance rejection under high wind loads. Since the key 
issue for implementing disturbance observer-based control is to 
generate proper compensation control signal when the types of 
the control signal and the disturbance signal are different, 
an optimal control model with a linear quadratic regulator 
(i.e., LQR-DOB) is constructed to precisely generate the 
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Table 5 
RMSs of the Position Errors Of LCT Antenna Using the Proposed DOB-CPC 
(H..+LQR-DOB) on Different Models Under Different Wind Disturbances 
(Time Duration: 10-20 s) 


Collaborative 
Average wind speeds and fluctuation Mathematical simulation 
amplitudes model testbed 
vy = 10 m s™!, Aan, = 2 ms! 0”3658 0”5081 
vw = 15 m s™!, Aav, = 6 m s™! 1”9860 40290 
vw = 20 m s™!, Aav, = 10 m s™! 7”8560 18”3300 


compensation control signal, which can make the actuator 
output proper physical quantity relatively close to the estimated 
disturbance and consequently compensate the disturbance as 
much as possible. Finally, to maintain the rapidity of the 
transient process and acquire the disturbance rejection 
performance provided by the LQR-DOB feedforward con- 
troller and the H,, feedback controller at the same time, a 
controller switching policy is designed to integrate the LQR- 
DOB feedforward controller, the H,, feedback controller, and 
the existing proportional controller as an observer-based 
composite position controller (DOB-CPC). 

To further verify the performance of the proposed DOB-CPC 
(H,.+LQR-DOB), we conduct simulation experiments, which 
includes mathematical simulation experiment and collaborative 
simulation experiment, for LCT’s pointing control system. The 
results of the mathematical simulation experiment indicate that 
the proposed DOB-CPC (H,.+LQR-DOB) performs much 
better than other controllers for all types of wind disturbances. 
Meanwhile, the results of numerical experiments on the 
collaborative simulation testbed (i.e., the collaborative simula- 
tion experiment), which integrates the multibody dynamics 
model of antenna constructed by Adams, further ensures that 
proposed DOB-CPC (H,.+LQR-DOB) is applicable in 
practice. 

Once the relocation and refurbishment of LCT is completed 
in 2024, the proposed DOB-CPC (H,.+LQR-DOB) will be 
verified physically during the astronomical observation to 
further evaluate its capability of achieving high precision 
pointing accuracy under large wind load. Based on this, our 
future work will focus on enhancing the integration between 
the proposed pointing control method (i.e., DOB-CPC) and the 
pointing correction technique to further improve the pointing 
accuracy of LCT. 
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